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Spin polarizabilities provide information on the internal structure of hadrons in the presence of 
weak external electromagnetic fields, and are actively studied by Compton scattering experiments. 
They provide finer detail than the regular polarizabilities since they require space and time-varying 
fields. Using an effective action in the weak field limit, we have identified methods to isolate 
each of the physical quantities (ji , a , j3 , Je\ , Jm\ , Yei , Jmi) for spin- 1/2 hadrons, both neutral and 
charged. We also perform a lattice QCD simulation to investigate the feasibility of the effective 
action approach. 
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1. Introduction 

The interaction of a spin- 1/2 hadron with a weak external electromagnetic field can be de- 
scribed by an effective, non-relativistic quantum mechanical (QM) action in Euclidean space as [1] 

Sqm = J d 4 xL QM , (1.1) 
where the effective QM Lagrangian is given by 

L QM = V f M (^^+{-N-qA) 2 -pio-B+^aE 2 -l L 2 

yf(x,t). (1.2) 



+ l -Y E \0-ExE- l -y M id-BxB- l -YE2 OiEijBj - y M2 OiBijEj 



Here 



E = ~-VA 4 , B = VxA, S = ^, E u = ^{ViEj + VjEi), etc. (1.3) 
and the various physical quantities are 

jj. magnetic moment 

a electric dipole polarizability 

/3 magnetic dipole polarizability 

Jei electric dipole to electric dipole (El— s-El) spin polarizability (1.4) 

Ym i magnetic dipole to magnetic dipole (Ml— >M1) spin polarizability 

Je2 magnetic dipole to electric quadrupole (Ml— »E2) spin polarizability 

Ymi electric dipole to magnetic quadrupole (El— >M2) spin polarizability 

This effective description of the interaction is expected to be valid for fairly weak fields. The 
polarizabilities encode rich information on the internal structure of the hadron with varying detail. 
The same polarizabilities appear in low-energy Compton scattering amplitudes [2], and are actively 
pursued by experiments and phenomenological studies. Although /x, a and j8 have been studied in 
lattice QCD [3], work on spin polarizabilities is just starting [4]. Our goal is to study the feasibility 
of isolating them on the lattice, combining effective theory and lattice QCD. 

2. Lattice discretization 

The two-point correlation function in the effective theory (or effective hadron propagator) 

G ss i(t,p,A u ) = axe' ; . (2.1) 

Since the effective action is in the bilinear form Sqm = Y^Ky, the path integration is easily per- 
formed and the resulting correlator G ss / is the inverse of the matrix K. After projection to finite 
momentum, the correlator is simply a 2 by 2 complex matrix in spin space, with s = 1,2 denoting 
spin up and down. 
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We discretize the effective action on a finite lattice of extent N x N y N z N t and spacing a, and 
evaluate the inverse numerically using a double-precision BiCGSTAB solver (since K is not her- 
mitian) with a convergence criterion of 10 -15 . Since the system is non-relativistic, only forward 
propagation is considered. The derivatives are replaced by appropriate differences on the lattice. 
We use the unit system in which h = c = 1 and e 2 = 1/137, and measure all other quantities in 
terms of fm. So the coordinates (x,y,z,t) are in fm, the nucleon mass m = 938/197 = 4.76 fin , 
the nuclear magneton jUjv = e/(2m) in fm, E and B fields in fm~ 2 , a and j8 in fm 3 , and y 's in fm 4 , 
etc. 

Most of the results below are obtained on a 24 3 x 64 lattice at a = 0.1 fm. The source location 
is at (12, 12, 12,4). We use Dirichlet boundary condition in the time direction, and study different 
boundary conditions in the spatial directions. This is one advantage of the effective approach: we 
can study boundary conditions exactly as they are in lattice QCD. Another advantage is that we can 
dial individual terms in the effective action and see their effects in the correlator. No such freedom 
is afforded in lattice QCD where all relevant terms are present simultaneously. Our basic strategy 
to determine the polarizabilities is: given a lattice QCD correlator, an effective QM correlator is 
computed by adjusting the polarizabilities until a match is found. 

Below we test the effective QM action with some input values in Table 1 (see [2]). We want to 
see if the input values could be recovered as output from the effective correlator alone. Such testing 
is useful because the same methodology can be used on lattice QCD conelators. The advantage 
here is that the effective QM correlator is fast to compute and free from Monte Carlo noise. Once 
a method is identified to isolate them, we can determine their true values by matching the effective 
QM correlator with the corresponding lattice QCD one. 



Table 1: Proton and neutron values for testing purposes. 



m 




a 


/3 


7ei 


Ym i 


7E2 


YM2 


(fm- 1 ) 


Ow) 


(l(T 4 fm 3 ) 


(l(T 4 fm 3 ) 


(l(T 4 fm 4 ) 


(l(r 4 fm 4 ) 


(l(T 4 fm 4 ) 


(l(T 4 fm 4 ) 


p 4.76 


2.79 


12 


2 


-3.4 


2.7 


1.9 


0.3 


n 4.76 


-1.91 


12 


2 


-5.6 


3.8 


2.9 


-0.7 



3. Free field 

In the case of periodic spatial boundary conditions, which we denote as bc=(l, 1,1,0) (the 
boundary condition in time is always Dirichlet in this study), the correlator projected to zero mo- 
mentum is a single exponential G(t) = wexp(— Xt) where the spectral weight is w = 1/(1 +am)/a 
and the mass of the particle is related to X by X = ln(l +am). This was confirmed numerically. In 
the case of Dirichlet in x and periodic in y and z, denoted as bc=(0, 1,1,0), the correlator contains a 
tower of states with discrete energies 

N, 2 

G(t) = Y w n e-^, where X n = ln(l +E n ) and E n = m+^-. (3.1) 

Here p n is the discrete lattice momentum given by 

sm(nna/(L x + a)) 

Pn = — -jz with n = 1,3,5, ••• , (3.2) 

a/2 
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where L x is the lattice size in x direction. Note that the even terms are absent under the specific 
boundary conditions. Since both w n and X n can be calculated exactly in the free-field case, we are 
able to confirm numerically that the correlator is indeed described by Eq. (3.1) and Eq. (3.2) to 
double precision. Similar checks were performed for other boundary conditions used in this work: 
bc=(0,0,l,0), bc=(0,l,0,0), and bc=(l,0,0,0). 

4. Constant electric field: a 

Next we turn on the electric field. The vector potential A = {—iEt ,0,0) in Euclidean space 
(we set A4 = 0) corresponds to a real-valued electric field E in the x direction in Minkowski 
space. We impose Dirichlet boundary condition in the x direction to eliminate the unphysical 
wrapping-around-the-lattice effects [5] present in lattice QCD so our boundary condition in this 
case is bc=(0, 1,1,0). For neutron, the interaction part of the effective action (all terms except the 
first two in Eq. (1.2)) is L\ nt = —\/2aE. We checked that a can be recovered straightforwardly 
from the correlator for different fields. For proton, the electric field causes an acceleration term 
in the correlator, in addition to the a term. We checked that a can be cleanly extracted from the 
unpolarized ratio G(t,a)/G(t,a = 0), which not only cancels out the acceleration term, but also 
the common energy term m + p 2 /{2m). The advantage here is that we do not need to know the 
details of the acceleration: its effects are fully incorporated in the effective QM correlator. Since 
the acceleration is the same in the effective theory and QCD, its effect can be canceled out in the 
same way in QCD in the ratio Gqcd((, g:)/Gqm(A 05 = 0). This is an alternative method to the one 
used in previous studies of charged hardons [6]. 

5. Constant magnetic field: /1 and /3 

Now we turn on the magnetic field. The vector potential A = (0,Bx,0) corresponds to a real- 
valued magnetic field B in the z direction in Minkowski space. The quantization condition ea 2 B = 
n(2n/N x ) would ensure constant B on the periodic lattice, but the fields produced are too strong: 
the lowest quantized B=306 fm~ 2 on our lattice would cause a proton energy shift of 1.5 GeV 
from piB (out of 0.938 GeV). One could reduce the strength from another quantization condition 
ea 2 B = n(2n/N x N y ) by 'patching' up the field on the edge of the lattice [7] to ensure uniform flux 
through every plaquette, but the lowest field (12.8 fm~ 2 ) is still too strong for our puiposes. So we 
abandon quantization and use weak fields (as small as B=0.2 fm~ 2 ). We impose Dirichlet be in x 
and put the source in the center of lattice to minimize the effects from breaking the quantization 
condition, and Dirichlet be in y to eliminate the unphysical wrapping-around-the-lattice effects, so 
bc=(0,0,l,0). For neutron, the interaction in the effective action is L in/ = —piBo z — 1/2/3 B 2 . We 
checked that pi can be recovered from the ratio of polarized correlators G\\(t) / G2i(t) and /3 can 
be recovered from the product G\\(t) * G2i{t). For proton, there is the additional complication of 
Landau levels. We checked that the extraction of pi is unaffected since the Landau-level effects 
cancel out in the ratio. For /3, the Landau-level effects can be cleanly removed by taking the double 
ratio 

Rp = [Gu(t,p)G 2 2(t,P)]/[Gu(t,P=0)G22(t,p=0)). (5.1) 



4 



Spin Polarizabilities on the Lattice 



Frank X. Lee 



6. Time and space varying fields: y 's 

For Ye 1, the choice of the vector potential 

A = (e\t 2 /a,ie 2 t,0), (6.1) 

where e\ and e 2 are real parameters, gives a time-varying electric field 

E = (-2eit/a,-ie2,0). (6.2) 

Because the vector potential has components in both x and y directions, Dirichlet boundary con- 
ditions are imposed in these directions to eliminate any wrapping-around-the-lattice effects. The 
boundary conditions are denoted as bc=(0,0,l,0). The interaction Lagrangian has the form 

L ira = \/2a{4t 2 e\/a 2 - e\) + y E \e\e 2 o z /a. (6.3) 

Due to the o z dependence Jei can be isolated from the ratio of polarized correlators 

R tki = Gu/G 2 2 ~ exp(2y El eie 2 t/a), (6.4) 

despite the complicated a term. For proton, the electric field causes an additional acceleration 
term, but the ratio is not affected and Jei can be extracted the same way. In both cases (neutron 
and proton), the method works out nicely numerically. We also checked analytic continuation: 
A = (e x t 2 /a,ie 2 t,0) and exg(2jE\_e\e 2 t / a) give identical results for Jei asA = (e\t 2 /a,e 2 t,0) and 
e\p(2iy E ieie 2 t /a) for small e\ and e 2 values. 
For Ye 2 , the choices are 

A = (e iy /2,0,ie 2 tz/a), E = (0,0, -ize 2 / a), B = (0,0,-«i/2). (6.5) 

Dirichlet boundary conditions are imposed in the x and z directions in this case so bc=(0, 1,0,0). 
The interaction Lagrangian has the form 

L int = \l2\xe y o z - l/2az 2 e 2 2 /a 2 - l/8j8cf + \/Ay E2 eie 2 a z /a. (6.6) 

The unpolarized a and j3 terms can be eliminated in the ratio G\ \ /G 22 which is left with both fx 
and Yei terms. Fortunately, the ji term is linear in the fields and Je 2 term quadratic. So the jj, 
contribution can be eliminated by averaging the mass shifts over (e\,e 2 ) and (—e\,—e 2 ). At the 
correlator level, the eliminations can be achieved simultaneously by the double ratio 

Ry E1 = [Gn{e\,e 2 )G\\(-e\,-e 2 )]/[G 22 (e\,e 2 )G 22 (-e\,-e 2 )] ~ &y^(\/2y E2 e\e 2 t /a). (6.7) 

The method applies to both charged and uncharged hadrons. It also implies that in order to deter- 
mine Je2 hi QCD we need to perform two lattice QCD calculations: one with an original set of 
fields, the other with the fields reversed. 
For 7M2, the field choices are 

A = (0,eit 2 /a,ie 2 t/2), E = (0,0, -ie 2 /2), B= {-e x xj a,0,e x zj a) (6.8) 
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Dirichlet boundary conditions are imposed in the y and z directions in this case so bc=( 1,0,0,0). 
The interaction Lagrangian 

Li„t = liei(xa x -zc z )/a- l/Saej- \/2fie\{x 2 + z 2 )/a 2 - l/4y M2 <?ie 2 tf z / a (6.9) 

has a similar form as that for Ye2, so the same methodology can be used to extract Jm2- Namely, 
the double-double ratio in Eq. (6.7) would be proportional to exp(— 1/27^2^1^2^/0)- 

Finally, for Jmu no choice of the vector potential can isolate it independently. We found that 
the choices 

A = (0, ie\tz/a,e2x), E = (0, —ie\z/ a,0), B = (—ie\t/a, —e2,0), (6.10) 
result in an interaction 

Li nt = pL{ieito x /a+e20 } )-\/2az 2 e\-\/2p(-e\t 2 /a 2 +el)-\/4(2yMi-yE2)e\e 2 o z /a, (6.11) 

that has contributions from both Jmi an d Je2- We can isolate the combination (2jmi — 7e2) in the 
same way as Jm2, then use the previously determined Je2 on the same lattice to pin down Ym\- 
Dirichlet boundary conditions are imposed in the y and z directions in this case so bc=( 1,0,0,0). 

In summary, these numerical tests demonstrate that all four spin polarizabilities can be disen- 
tangled by a judicious choice of the vector potential, and manipulation of effective correlators. The 
real challenge is whether the same can be done with lattice QCD correlators which are born with 
Monte Carlo noise. 

7. y E \: a lattice QCD case study 

Here we perform a preliminary real lattice QCD simulation to see whether the strategy pro- 
posed is feasible. The QCD data are generated on 24 3 x 48 lattice with Wilson actions at a=0.093 
fm. The field is applied by A = (e\t 2 /a,e2t,0) and bc=(0,0,l,0), with field values (^1,^2) = 
(0.23,0.46) fm~ 2 . We analyzed 700 configurations for 6 pion masses from 893 to 404 MeV. In 
Figure 1, we see that the nucleon ground state starts to dominate at t = 7. The effective mass of 
Re [Gn (t)] is very different in neutron and proton (mostly due to acceleration effects of the proton), 
but the ratio 1m[G\\{t) / G22{t)\ ~ sm[yE\e\e2t /a] is the same. The same ratio in QCD suggests a 
negative Je\ for both neutron and proton, but suffers from large statistical uncertainties. The 3 
lines are predictions at 3 different values of Jei in units of 10~ 4 fm 4 . The QCD data suggest a wide 
range for Je\ (anywhere between and -200). To get a more definitive value for y^i , the error bars 
need to be reduced by at least a factor 10, which means the statistics need to be improved by 100 
times if no other improvements are made. 

8. Conclusion 

The effective QM correlator provides a controlled and fast way of examining spin polarizabil- 
ities and their systematics on the lattice. By designing the background fields and manipulating the 
correlators, we have identified methods to isolate all the physical quantities {\i , a, j8 , Yei , Ymi , 7e2, 7m2) 
for spin- 1/2 hadrons as defined in Eq. (1.2). This is possible because of the different field and spin 
dependences in the terms containing these quantities. By matching with the corresponding lattice 
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Figure 1: Some results from the Je\ study. 

QCD correlators on the same lattice with the same boundary conditions, their QCD values can 
be determined in principle. A preliminary study of y%\ using lattice QCD data suggests that the 
challenge lies in the Monte Carlo noise present in the lattice QCD simulations: the noise must 
be reduced significantly before reliable information can be extracted. Several strategies are being 
explored to this end. 
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